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Abstract 

This paper aims to cast some new light on controlling chaos using the OGY- and 
' the Zero-Spectral-Radius methods. In deriving those methods we use a generalized 

. procedure differing from the usual ones. This procedure allows us to conveniently 

treat maps to be controlled bringing the orbit to both various saddles and to sources 
with both real and complex eigenvalues. We demonstrate the procedure and the 
subsequent control on a variety of maps. We evaluate the control by examining the 
basins of attraction of the relevant controlled systems graphically and in some cases 
analytically. 

> 

q ! 1 Introduction 

One of the basic properties of a chaotic dynamic system is topological transitivity. We 
can take advantage of that property. In practice it can be used to get from any point 
in the attractor to the vicinity of a given point. The popular OGY-method is typically 
used when the target point is an unstable fixed point or a point in an unstable periodic 
cycle. Once we have arrived close enough to the target point we use linear control with 
respect to a parameter in order to stay in that neighborhood. The OGY-method and 
its various versions aim at bringing the orbit to the stable set of the target point and 
keeping it there. In that way the system itself will attract the trajectory to the target. 



C3 



The method is described in [14] and among many other well formulated descriptions and 



derivations we want to mention 0, ||, Among others, |l5| deals with special techniques 
for stabilising unstable periodic cycles of higher order. Several other variants of local 
control through parameter perturbation and their comparison to the OGY-method are 
presented in [HJ]. An investigation of the relations between these methods and the 

"pole placement" methods can be found for instance in [|13[. An attempt to establish a 



general methodology is made in || 11, [12]. In order to get close to a certain point in a 
chaotic attractor rapidly from a relatively distant point,various versions of the "Targeting"- 
method can be used. The original idea of that method is described in [[!(], 17, 18]. A large 
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overview of early chaos control methods was made in And recently the same authors 
have assembled most of recent work on the subject in ||. 

In this paper methods of controlling by parameter perturbation in a map are used to 
stabilize the system to an unstable fixed point. A complementary set of examples can 
be seen in [3]. We will use the OGY-method and the Zero-Spectral-Radius method. The 
examples were chosen in order to get a broad variation in the properties of the unstable 
fixed points or periodic cycles. They vary from saddles to sources. Some eigenvalues of 
the corresponding Jacobian matrices are real and some are complex. We investigate two- 
and three-dimensional maps. 

In some cases control is applied several times on the same map, each time perturbing 
a different parameter. 

Our main goal is to evaluate the efficiency of the above methods on the variety of 
maps they are tested. The efficiency can be viewed through the plots (with level curves) 
of the basins of attraction of the controlled systems. In the case of the three-dimensional 
maps the basins of attraction are viewed in perpendicular cross-sections through the fixed 
point. Explicit expressions for the controlled systems are derived in several cases. In some 
of the cases the basin of attraction of the controlled system is so small or so narrow that 
the local control method in question is to be employed with special precaution or to be 
preceded by a "targeting" procedure. 

This paper is organized in the following manner: 

Section |2|. The derivation of formulae for linear control by parameter perturbation. 
Both the Zero-Spectral-Radius method (ZSR) and the OGY-method are derived through 
matrix manipulations. 

Section |3|. A study of the control of a 2-dimensional discrete map: 

The Henon map is stabilised to a saddle fixed point in several cases. Control is imple- 
mented with respect to several available parameters and using the two methods mentioned. 

Section ^. A study of the control of three-dimensional discrete maps: 



In subsection 4.1 a Henon- like map is controlled to a fixed point with one unstable real 



eigenvalue and two stable complex ones. Both the OGY- and the ZSR-methods are used. 



In subsection [4.2| a system of three coupled logistic maps is controlled with the OGY- 



method. The map is stabilised to a fixed point with one stable and two unstable directions. 

2 General theory and control formulae 

Let us consider the map 

x^F(x,p), (2.1) 

where x is an n-dimensional vector and p is a real parameter. Without loss of generality 
we assume the fixed point x* to be the origin for the parameter value p = (if we are 
interested in other periodic orbits we here use the corresponding iterate of / instead of / 
itself). The map can then be written in the form 



x^Jx + pw + Fj, (2.2 

e J is the Jacobian matrix and w = ^ , 
origin) and with p = 0. Fj is of higher order. 



where J is the Jacobian matrix and w = both calculated in the fixed point (here: the 
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To control the system we choose an expression for p being a linear function of the 
components of x : p = ofjX. 

Thus the function representing the controlled system can also be written in the form 

x^J c x + Fj, (2.3) 

where J c = J + wofj. 

The next step is a change of coordinates. Suppose the new coordinates u are defined by 



x = Qu, where Q is an n x n non-degenerate matrix. In the new coordinates the map (2.1) 
can be written in the following form: 

u -» Mu + pb + F M , (2.4) 

where M = Q~ l JQ, b = Q~ l w and Fm is of higher order. 
The controlled system can be written in the form 

«-»• M c u + Fj, (2.5) 

where M c = M + baj^, = ctjQ. 

The case where M becomes a diagonal matrix or a matrix in Jordan normal form is of 
special interest. 

We shall look at ways of choosing the coefficients in order to make the control work. 
If the coefficients in Um (equivalently aj) are chosen so that the absolute value of the 
eigenvalues are less than one the fixed point becomes stable for the controlled system. 
This means that we can control the system in the vicinity of the fixed point (actually the 
controlled system will converge towards the fixed point). 

We shall use two methods: the Zero-Spectral-Radius method (ZSR) and the OGY- 
method. In the ZSR-method all eigenvalues of the Jacobian matrix for the controlled 
system are required to be zero. In the OGY-method only the unstable eigenvalues, that 
is the ones with absolute value greater than one, are required to be zero and the stable 
ones, those with absolute value less than one are left unchanged. Let's have a look at the 
first case. 



2.1 The Zero-Spectral-Radius method (ZSR) 

We can derive general expressions for the coefficients on (the components of %) making 
all the eigenvalues of M c equal to zero. 

Theorem 2.1. Let the matrix M be a diagonal matrix 



( A1 


• 


• 


\ 





A 2 • 


• 




V o 


• 


A n 


/ 



where Aj ^ A,- for i ^ j. Furthermore we suppose that all bi ^ 0. Then the eigenvalues of 
M c are all equal to zero if and only if 

—\ n 

Ui = — *- , i = l,...,n. (2.6) 

n (A? — Xj)bi 

(Aj can be complex). 
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Proof. The characteristic polynomial P(X) = det(M c — XI) = (— l) n A n , since all its roots 
are equal to zero. We will develop M c = M + ba M . We get 



Mr 



( Ai + b\oc\ b\a 2 
b 2 a\ A2 + b 2 a 2 



\ b n ai 



b n a 2 



b\a n 
b 2 a n 

An + b n a n j 



and P(X) = det(M c - XI) equals 



det 



/ Ai + b\ui 
b 2 ct\ 



X 



\ 



b n ai 



b\a 2 
X 2 + b 2 a 2 

b n a 2 



biot n 
b 2 a n 

X n + b n a n - X j 



Now develop P{X\) by elementary column operations, which results in 



ot\ det 



b 2 



b\a 2 
X 2 — X\ + b 2 a 2 



b n a 2 



a\ det 



/ 61 

b 2 

\ b n 

n-1 





A2 — Ai 




b\a n 
b 2 a n 



\ 



An - Ai + b n a n j 



\ 



An — Ai / 



(Ai - A 2 )(Ai - A 3 ) • • • (Ai - A„) = {-\) n X n x . 

1. For ai, i = l,...,n we get the analogous 



= b\ai{— 1 

Solving for a\ we get (|2.6|) for i 
expressions as necessary conditions. 

If the characteristic polynomial has coefficients Cj, i = 0, . . . , n, that is P(X) = c n X n + 
• • -+ciA+co, then we can immediately see from the determinant expansion that the leading 
coefficient must be (— l) n then giving a system of n equations c n -iA" _1 + - • -+ciAj+co = 0, 
i = 1, ... ,n which has only the trivial solution because the Aj:s are all supposed to be 
different from each other. (The system can be written in the form C(c n -i, • • • , c\, cq) t = 
and it is well known that in this case the determinant of C is different from zero). 
Thus the conditions are also sufficient. The proof is finished. 

We shall now look for the concrete formulae in some special cases used in this paper. 

Suppose J is two-dimensional and has two real eigenvalues Ai 7^ A2. Then we get the 
following formulae for the coefficients: 



a 1 



-A? 



6i(Ai-A 2 )' 



OL 2 



\ 2 
~ A 2 



b 2 (A 2 - Ai) 



(2.7) 



Suppose J is two-dimensional and has two complex eigenvalues A ± i/i, where [i 7^ 0. 
We remark that 



S 



-1 



A + fii 






X — fj,i 



S 



X ^ 
—fj, X 
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where 

S = 



i 1 
1 i 



Thus the expressions for a.\ and a 2 in this case can be obtained from the case with two 
different real eigenvalues if we put Ai = A + X 2 = A — \xi and use b = Sb and a T = a T S, 
where b gives b\ and b 2 for the real eigenvalue case and b gives b\ and 62 for the complex 
eigenvalue case and analogously for a and a. 

Consequently if M is of the form 

A n 
—fi A 

where ^^0 then we get the following expressions for the coefficients: 

«i = Mfi 2 - A 2 ) - 2b lf iX)/A, a 2 = MX 2 - fi 2 ) - 2b 2 »\)/A, (2.8) 

where A = fi(bf + b 2 ). 

Suppose now J is three-dimensional and has one real eigenvalue A u and two complex 
eigenvalues A ± ifi, where /j^O. Suppose M is of the form 




where Then we get the coefficients: 

ai = -Xl/iDh), a 2 = -(Bb 2 + Cb 3 )AD, a 3 = (Cb 2 - Bb 3 )/AD, (2.9) 

where 

A = n{b\ + bj), B = 2X 3 fi + X u fi 3 + 2^i 3 A - 3X u X 2 fi, 

C = A 4 - A U A 3 + 3X u fi 2 X -n 4 , D = X 2 + /j 2 — 2X U X + X 2 U . 

2.2 The OGY formulae 

The same techniques and the formulae in the previous subsection can also be used for the 
OGY-control. We give the result in some cases. 

Suppose for a two-dimensional system the fixed point is a saddle and M has the form 

Ai 
A 2 

where the absolute value of Ai is less than one and the absolute value of A2 is greater than 
one. Then we get the coefficients: 

ai =0, a 2 = -A2/62. (2.10) 

This is easily seen to coincide with the known form of the OGY- formula. We have 
a J = a^Q^ 1 . And since the columns of Q consist of the eigenvectors of J the rows of 
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rp 

Q^ 1 consist of the contravariant basis vectors f\ and f 2 . This leads to aj = — A2/2 /&2- 
But because b = Q~ l7 w we get 62 = f % ■ resulting in the familiar 

— T 

_T ^2/2 

aj = 

Suppose for a three-dimensional system the fixed point is a saddle with two-dimensional 
unstable manifold and M has the form 




where the absolute values of Ai and A2 are greater than one and the absolute value of A3 
is less than one. 

Then we get the coefficients: 

aX = W^MY a2 = RA^)' 03 = (2 ' U) 
similar to formula vtl. 



3 Comparing ways of controlling a two-dimensional map 

We have used our method to control some two-dimensional maps with only linear and 
quadratic terms. The fixed points were saddles and even sources. In this section, however, 
we will concentrate on showing different ways of controlling the familiar Henon map with a 
saddle. The cases with sources are investigated in Q and shortly treated in the discussion. 



3.1 The Henon map 

Let us now consider the map 

J\ / a-bY -cX 2 



YJ V dX + e ' ,:U » 

For c = d= l,e = this is the famous Henon map. The parameters c, d and e are 
introduced in order to apply control in several ways to the standard Henon map. Control 
is applied by perturbing one of the parameters while the other ones are held fixed. For 
instance applying control with respect to the parameter a is carried out by permitting the 
parameter to vary about a nominal value a*, so that a = a* + p. Applying control with 
respect to any of the parameters b, c, d and e is carried out analogously by adding the 
perturbation p to the nominal value of the parameter in turn to be used. In addition we 
apply control by adding multiples of p to both coordinates. 

In the following, c* = d* = 1 and e* = will be used. We will assume that the values 
of a* and b* are such that the map has a chaotic attractor and that within the attractor 
there is a saddle fixed point (x*,y*), with one stable direction with the stable eigenvalue X s 
and one unstable direction with the unstable eigenvalue X u . Notice that, given the above 
values of the parameters, x* = y* . 
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One example of such a case is the common version of the Henon map with a* = 1.4 
and b* = —0.3. It is treated in many papers and textbooks. In the numerical example in 
Section |3.1.7 we will consider a special case with values a* = 1.05 and b* = —0.5. Iterating 



the map yields a familiar multi- folded horseshoe-shaped attractor shown in Fig. ||a). One 
reason for choosing this particular version of the Henon map is that the stable eigenvalue 
of the fixed point is larger (0.2665) than in the "common" case (0.1559). This fact leads to 
bigger differences in the results when we apply the two control methods, ZSR and OGY. 

We control the system by the ZSR and OGY methods. In order to use the formulae in 
Section § we transform the system in two steps: 

1) Moving the fixed point to the origin. Separating the linear and non- linear parts. 

2) Taking the eigenvectors as new basis vectors. 

3.1.1 Moving the fixed point and separating matrix terms 

We introduce the new variables (x,y) = (X — x*,Y — y*). The map (^) can now be 
rewritten as 

x \ I a — b(y + y*) — c(x + x*) 2 - x* 
y J V d(x + x*) + e - y* 

We apply control with respect to each of the parameters a, b, c, d and e one by one 



(3.2) 



as described in the beginning of Section [3J. Expanding and separating the matrix terms 
(keeping in mind that x* = y* , c* = 1, d* = 1 and e* = 0), we obtain the formula ( |2.2| ) in 
each case with the following expression for J: 

, , -2x* -b* 



1 

The expressions for w will vary: 

a)W=(l), b)w=(~f I. (3.3) 



v / V 

c)w= ( " ( ^ )2 y d)w=(^ or e)w=(j), (3.4) 

depending on which parameter to be perturbed. 

In the following we consider especially the first two cases a) and b). On the other hand 
we have applied control perturbing both a and e at the same time in the following cases 
(meaning (a, e) — > (a, e) T + pw, and control with respect to p): 

f) w = ( \ \ 9)w=(\\ and h)w=(^\ (3.5) 

Out of these, we especially consider the case f). 

The expressions for Fj is the same in cases a), e) , f) , g) and h) : 

The expression for Fj in case b) is 
Tj= f- P y-x* 
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3.1.2 Changing basis vectors 



The eigenvalues of J are A = —x* ± y (x*) 2 — b* which results in A1A2 = b* and Ai + A2 = 
— 2x* . As mentioned in the beginning of this section we are particularly interested in 
the fixed point being a saddle. The eigenvalues can then be noted X s and X u and the 
eigenvectors are 

A s \ , — ( X u 

1 I , and e u = I 

The basis is changed so that the eigenvectors of the Jacobian matrix become the new 
basis vectors. The coordinates in the new eigenvector basis are called (u,v), with u being 
the coordinate in the stable direction and v the coordinate in the unstable direction. The 
matrix Q is the transition matrix 

_ / A s A n 
V 1 1 

leading to a diagonal matrix M = Q~ l JQ. Since u = Q~ l x, we get 

(x-X u y), v = - 1 (-x + \ s y). (3.6) 



A s - A„ A. s — A 



Furthermore, since 

I ) - ^ 

we get for control with respect to a 

^s ~ A u A s — A u 

we get for control with respect to b 

b-> (3.8) 



As — A u As — X u 

and for control in the /)-case: 

= 1 + A M _ -1-As 

As — A u As — X u 



(3.9) 



3.1.3 Using Zero-Spectral-Radius control 

Although the fixed point in our Henon attractor typically is a saddle point it does not 
matter in the case of the ZSR-method if there is a stable eigenvalue or not. So we will here 
use the more general notation (Ai, A2) instead of (X s , X u ). We can now obtain expressions 
for a± and 02 in accordance with fl2.7p leading to explicit expressions for the parameter 
perturbation p = a±u + a<iv for the Henon map (using A1A2 = b and Ai + A2 = —2x*) 
In the case of control with respect to a we get 



Qi = —Xf, «2 = — A| and p = 2x*x + b*y. 



(3.10) 
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In the case of control with respect to b we get 



X 2 A| b* 
a i = ~T' °2 = — and p = -2x -y. 



(3.11) 



In the case f) we get 
-A? 



a i 



1 + A 2 



«2 



\ 2 
~ A 2 



1 + Ai 



(3.12) 



In the case f) the analytical expression for p cannot be brought to a simple form, so it 
is left out. Thus in the following only the a)- and 5)-cases are treated analytically. 



3.1.4 Deriving the controlled system for the ZSR-method 

Control parameters are determined from the linear part of the system, but nevertheless the 
non-linear part is still present and together with the parameter perturbation a new system, 
the "controlled system" is created. In general, the controlled system will differ more from 
the original one the further one gets from the fixed point. We insert the expressions for p 
into equation ( |2.2D and get the following controlled systems. 



In the case of control with respect to the parameter o: 



-x 

X 



(3.13) 



In the case of control with respect to the parameter b: 



( 2 ( A l + A2) A1A2 2 

x : zyH — y 



y 



y 



(3.14) 



which is the same as 



/ 2 b* o 

' -x 2 + 2xy + — y 2 

y 
\ 



(3.15) 



In the first case one can clearly see that the controlled system is only dependent on x 
and convergence to the fixed point (the origin) will occur under the same circumstances 
as for the one-dimensional map x — > — x 2 . So there is convergence for ( 3.13 ) to the fixed 



point (0, 0) 1 , that is (X, Y) 1 = (x*,y*) 2 , from starting points with X e]x* — l,x* + 1[. 
For X = x* ± 1 there is convergence to the point (x* — 1, x* — 1) T and divergence otherwise. 
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3.1.5 Using OGY-control 



Since the OGY-method is based on the very property of a fixed point being a saddle with 
at least one stable direction we will here use the notations (A s , X u ) instead of (Ax, A2). 
According to ( [2,10 ) we get p = a\u + a<iv = Ou 



fey. Thus 



V 



-X u v 



(3.16) 



Using (3.6) and ( |3.7| )~(^S|) we can get the particular expressions for each version of 
parameter perturbation. 



3.1.6 Deriving the controlled system for the OGY-control 



Inserting the expressions for p into equation ( |2.2D we get the following expressions for the 
controlled system applying OGY-control. 



In the case of control with respect to a: 



X s x — x 2 
x 



(3.17) 



In the case of control with respect to b: 



I 



\ 



X q x — x 



Xu . X S X U 2 
—xy H — y 

y y 



(3.18) 



Here again in the first case the controlled system is only dependent on x and conver- 
gence to the fixed point (the origin) will occur under the same circumstances as for the 



XcX 



So there is convergence for ( 3.17] ) to the fixed point 



one-dimensional map x - 
(0,0) T , that is (X, Y) T = 

X = x* + X s — 1 and X = x* + 1 there is convergence to the point (x* + A s 
and divergence otherwise. 



(x* , y*) , from starting points with X e]x* + A s — 1, x* + 1[. For 

l,x* + X s -l) T 



3.1.7 A numerical example 

The Henon map with parameter values a = 1.05 and b = —0.5, (and c = d = 1, e = 0) will 
be used. The corresponding chaotic attractor is shown in Fig. |l]a). The unstable fixed 
point is (0.8048, 0.8048) T and it is marked with a circle (o) in the plot. The eigenvalues for 
the Jacobian matrix are A s = 0.2665 and A n = -1.8760. The corresponding eigenvectors 
are e s = (0.2665, 1) T and e u = (-1.8760, 1) T . The two methods, OGY and ZSR, yield 
almost congruent basins of attraction using the same parameter (parameters b, c, d or e). 
On the other hand depending on the parameter to be perturbed the basins of attraction 
are very different in shape. In the cases g) and h) the basins are similar in shape and 
geometry regardless of the method (mostly almost vertical contour lines). The /)-case 
application of the two methods yields a vastly differing inner geometry, as can be seen in 
the Fig. go) and |c). In all cases the ZSR-method displays larger regions of the fastest 
convergence than the OGY-method perturbing the same parameter. 
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> 




> 




Figure 1: a) OGY-control in the case /. b) OGY-control with respect to the parameter b. 
c) ZSR-control in the case /. d) OGY- and ZSR-control with respect to the parameter a. 

In Fig. a), b) and c) the map is iterated 5 times. The central dotted regions contain 
the points from which the map converges and reaches within a distance of less than 0.001 
from the fixed point (o). The peripheral dotted regions represent points from which the 
map certainly diverges. 

Fig. a) displays the result of OGY-control in the /)-case. The level curves are: 0.001, 
0.01 and 0.1. Fig. b) shows the result of OGY-control with respect to the parameter b. 
The level curves are: 0.003 and 0.06. Fig. c) displays the result of ZSR-control in the 
/(-case. The level curves are: 0.00001, 0.0001, 0.001 and 0.01. Finally in Fig. d) the 
results of control using both OGY- and ZSR-methods with respect to the parameter a are 
displayed. It was stated in subsections 3.1.4 and 3.1.6 that the controlled systems only 



depend on x. The graphs show the number of iterates needed to get within a distance of 
0.001 from the fixed point from various starting points. The solid graph represents the 
ZSR-method. The dotted lines X = x* db 1 mark the boundaries of the basin of attraction. 
The dotted graph represents the OGY-method. The dotted lines X = x* + X s — 1 and 
X = x* + 1 mark the boundaries of the basin of attraction. 

In this section we have not so far paid particular interest in limiting the size of the 
parameter perturbation p. We have looked for the entire basin of attraction regardless of 
how much an unlimited value for p might distort the system. This calls for a discussion. 
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Where is the region of the smallest initial p:s? Let us have a look at the OGY-method 
first. The parameter perturbation p is generally expressed in formula ( 3.16| ). Clearly the 



smallest values for p are required close to the straight line determined by v = (formula 
(|3.6| )), that is along the direction of the stable eigenvector (X S ,1) T or (1,3.75) T . This 
applies in all of our cases and is no big surprise giving the basic idea of the OGY-method. 
Investigating the expressions for p (formulae ( 3.1C|) - (|3.11 ) and the other cases c)~h) not 



explicitly given in this text) in the ZSR-method yield correspondingly straight lines as 
the centre lines of stripes constituting the regions of the smallest initial p:s. These lines 
have directions (1, —2x*/b*) T or (1,3.22) T in the cases a), b) and c). The cases d) and e) 
have a centre line with direction (1,3.84) T . In many cases the direction of the line fits 
well in the image of the basins of attraction. For instance, in Fig. [!]&) this stripe of small 
p:s follows narrow parts of fast convergence in the basin of attraction. Nevertheless, the 
further one gets along the line from the fixed point the more powerful the non-linear part 
of the controlled system becomes yielding uncertain results. 

4 Controlling three-dimensional maps 

4.1 Controlling a 3-dimensional Henon-like map 

Let us consider a three-dimensional map somewhat similar to the Henon map: 



(4.1) 



where normally p = 0. 

We will investigate the effect of controlling the map with the OGY-method and the 
Zero-Spectral-Radius method. We will consider a case where there is a chaotic attractor 
including an unstable fixed point (the origin) with one unstable direction, with the eigen- 
value A u , and a stable manifold with two complex eigenvalues A2 and A3. We will apply 
control using the parameter p. 

4.1.1 The control formulae 




In order to write the OGY-control formula we refer to Section g formula ( 2.10 ). Analo- 
gously with that formula we have a\ = —X u /bi, a.2 = and 03 = yielding;? = —\ u u\jb\, 
where u\ is the first coordinate of the new basis vector u. 

As to the ZSR-control formula we refer to Section [2] formula ( |2.9| ). The derivation of 
the coefficients is quite tedious but is conveniently carried out with the aid of a computer. 

4.1.2 The controlled systems 

The controlled system of the map applying the control methods can be derived as follows. 
Consider the Jacobian matrix of the controlled system: 
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where A = a + cti, B = b + a.%, C = c + 0:3 and p = a\x + a^y + 032:. The characteristic 
equation is 

A 3 - AX 2 - BX - C = 0. (4.2) 

In the case of OGY-control the characteristic equation can also be written in the fol- 
lowing form since the one unstable eigenvalue is required to be zero: A(A — A2) (A — A3) = 0. 
Developing and comparing the two equations yield: A = A2 + A3, B = — A2A3, and C = 0. 
So the map has a controlled system of the following form when A u is required to be zero 
in accordance with the OGY-method: 

(A 2 + X 3 )x - X 2 X 3 y - x 1 

|. (1.3) 

y 

This implies that the controlled system does not depend on z. 

In the case of ZSR-control the characteristic equation can also be written in the follow- 
ing form since all eigenvalues are required to be zero: A 3 = 0. 

Comparing with the equation |4.2| yields: A = B = C = 0. So the linear part of the 
first equation vanishes and the map has a controlled system of the following form: 

(4.4) 
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As for the Henon map the map 4.4 will depend only on x and convergence will occur 
under the same circumstances as the one-dimensional map x — > —x 2 . That is the map 
will converge to the fixed point (0, 0, 0) T from any point with x g] — 1, 1[. For points with 
x = ±1 there is convergence to the point (—1, —1, — l) T and divergence otherwise. 

4.1.3 A numerical example 

We will consider the particular case: 

-1.65x-0.3y-0.2z-x 2 +p 







■ 


H 







x I . (4.5) 

y 

The corresponding attractor can be seen in Fig. |2]a). 

As can easily be seen the map has a fixed point in the origin (0, 0, 0) situated in the 
chaotic attractor. The eigenvalues and eigenvectors of the Jacobian matrix in this point 
are 

X u = -1.5395, 

e u = (0.7906, -0.5135, 0.3336) T , 
A 2>3 = -0.0552 ± 0.3562?, 

e 2 3 = (-0.1211 T 0.0076i, -0.0306 ± 0.3352i - 0.9059 T 0.2265i) T . 

That is, there are one real unstable eigenvalue X u and two complex stable eigenvalues 
A2 and A3. 
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Figure 2: a) The attractor. b) A cross-section of the basin of attraction through the fixed 
point at z = using OGY-control. 

We shall now have a look at the basin of attraction and the rate of convergence of the 
control methods depending on the starting point. Three cross-sections were made through 
the fixed point: those were the planes x = 0, y = and z = 0. The controlled system was 
iterated five times starting from each point on each plane. Contour lines show how close 
the system point has reached the fixed point after five iterations. We choose to display 
the result only from the cross-section z = 0. 

Controlling with the OGY-method. Applying the method in order to control the 
map to the fixed point (0, 0, 0) T works well within certain limits of distance from that 
point. First, the planes x = and y = only display parallel contour lines which is 
explained by the controlled system not depending on z. The contour lines of the plane 
z = are shown in Fig. |2|6). 

The fixed point is marked with a cross (+). The vaguely boomerang shaped contour 
lines represent the levels (counted from inside out): 0.005, 0.010, 0.020, 0.050, and 1.000. 
The last level represents more or less the boundary of the basin of attraction for the 
OGY-controlled system. The 0.020 level makes an extra turn at the left side. 

Controlling with the Zero-Spectral- Radius method. The y = and z = planes 
display parallel lines and x = does not display any contour lines at all. This is explained 
by looking at the extremely simple expression for the controlled system. 
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From all starting points in the x = plane the system goes exactly to the origin in the 
second iteration at the latest. 









l-( 








4.2 Coupled logistic map with OGY-control 

In this subsection we examine the efficiency of the OGY-control for another three-dimen- 
sional map. The map is stabilized towards a saddle fixed point with two unstable eigendi- 
rections. Another version of this map with a saddle point with one unstable eigendirection 
is examined in Q . We consider the system consisting of three coupled logistic maps 

7*1 x (1 — x) 

r2 y (1 - y) 
r 3 z (1 - z) 

where ri, 7*2, ^3 and p are parameters. The coupling parameter p is chosen to be the control 
parameter. As the state space of the system we consider the unit cube [0, 1] x [0, 1] x [0, 1]. 
With parameter values 7*1 = 3.9, 7*2 = 3.95, r% = 1.0 and p = 0.0736, the point 

x* T = (0.7291,0.7323,0.2889) 

is a fixed point and there are one stable and two unstable eigendirections. 

The stable eigenvalue and the unstable ones with the corresponding eigenvectors are 

A s = 0.35600 e s = (0.021204, 0.020760, 1.3728) T 

\ ul = -1.4094 e ul = (0.63971, -0.54033, 0.0063062) T 

\ u2 = -1.6750 e u2 = (0.57127, 0.65842, 0.080606) T . 

The vector w becomes 
df 

w = -J-(x*,p*) = (-0.56084, -0.57308, 1.1339). 
op 



Here we calculate the control coefficients according to ( p. 11 ) using again in Section 2.2 
x — x* and p — p* = bp instead of x and p. 

The explicit expression for the control is obtained as 

bp = 104.33x - 107.46y + 0.013503z + 2.6172. (4.6) 

The controlled map has the Jacobian matrix 

/ -60.036 60.130 0.023502 
J c = -59.921 60.016 0.023337 
\ 118.17 -121.98 0.37538 

with (0.11687,0.1163,1) as eigendirection corresponding to the eigenvalue zero and 
(1,0.97904,64.737) as eigendirection for the non-zero eigenvalue. 

The cross-sections through the fixed points of the small region of convergence are de- 
picted in Fig. ||. The convergence levels are (calculated using 10 iterates) (0.048, 0.125, 0.5, 
1,5,10,50,100) • 10~ 9 . 
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0.7288 0.729 0.7292 0.7294 



0.732 0.7322 0.7324 0.7326 



Figure 3: One stable, two unstable eigendirections. Convergence levels for a) y = 0.7323 
and b) x = 0.7291. 



5 Discussion 

The message of this paper is two-fold: 

1 . The display of basins of attraction with level-curves for the controlled systems for a 
variety of maps. 

2. A new procedure for deriving control formulae for the OGY-method and related 
methods. 

We have chosen a variety of maps covering some unusual cases with real or complex 
eigenvalues. For this purpose we have brought forward a procedure, which in our opinion 
simplifies the derivations of control formulae through matrix diagonalising. Moreover the 
procedure yields convenient expressions (Section ||). The OGY-method and the ZSR- 
method thereby appear as mere special cases of a general way of thinking, regardless of 
the number of dimensions, the proportion of unstable eigendirections or the quality of 
the eigenvalues. There are possibilities to bring forward other similar methods than the 
mentioned ones. A comparison with the original form of the OGY formula is given in 



subsection |2.2| . Our procedure is implemented in all sections but it is displayed more in 
detail in Section ||. This section ^ contains a broad study of many possible ways to control 
the well-known Henon map. The goal of the control is to bring the orbit to a saddle fixed 
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point, which works successfully in all cases. But the geometry of the basin of attraction 
varies considerably depending on how the control is applied. 

Our "three-dimensional Henon" map has two complex eigenvalues in the origin and 
serves as an example of applying our procedure in such a case. 

In the case of the three-dimensional coupled logistic map all the eigenvalues are real, 
one of them stable and two unstable. The control coefficients become very large and the 
basin of attraction to the fixed point very small. The eigenvector corresponding to the 
eigenvalue zero (in the controlled system) dictates the direction of fastest convergence. 

Further examples with convergence level graphs can be found in Q: 

• Another version of the coupled logistic map with one unstable and two stable eigendi- 
rections and controlled with the OGY-method. 



Two-dimensional maps are stabilised to sources with either real or complex eigen- 
values. Obviously, the OGY-method is not applicable. The control with the ZSR- 
method worked well although the basins of attraction were rather narrow due to the 
fact that the two eigenvalues were near each other in size. This was the case also in 



Section 4.2. 



• There is also an example on OGY-contr oiling the Poincare map of a Duffing system 
either to a saddle fixed point or to a two-periodic saddle point. An analysis of the 
rate of convergence is included. 

As to the local speed of convergence: the ZSR is always the more efficient method due 
to the fact that the eigenvalues of all directions are placed at zero. The speed of the OGY 
method depends on the value of the stable eigenvalues. If their absolute value is close to 
one, the convergence is slow. If it is close to zero the behaviour is similar to the one of the 
ZSR-method. On the other hand, further away from the fixed point the importance of the 
non-linearity of the map grows. There the properties of the linearisation almost vanish. 
Therefore, in general, both methods display almost congruent outer boundaries of their 
basins of attraction. The non-linear influence also create fairly distant areas with fast 
convergence requiring only limited control. Much closer areas may have slow convergence 
or even divergence. 

We have concentrated on applying linear control to one parameter at a time. In general 
we would archive better results by perturbing several parameters or by applying non-linear 
control. Another procedure would be to combine the linear methods with targeting. 
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